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Abstract 

A  novel  three-dimensional  hybrid  Finite  Element  Method  (FEM)  and  Physical  Optics  (PO)  tech¬ 
nique  for  the  efficient  analysis  of  scattering  problems  is  presented.  It  makes  use  of  FEM  for  the 
regions  with  small  and  complex  features  and  PO  for  the  analysis  of  the  electrically  large  objects  of 
the  structure  taking  into  account  mutual  interactions  between  the  FEM  domains  and  the  objects  ana¬ 
lyzed  with  PO.  The  hybrid  method  proposed  makes  use  of  an  iterative  FEM  for  open  region  problems 
that  allows  the  FEM  domain  to  be  truncated  with  a  minimum  number  of  unknowns  while  retaining 
the  original  sparse  and  banded  structure  of  the  FEM  matrices,  and  also  an  easy  hybridization  with 
other  numerical  techniques  (as  PO  in  this  paper).  Several  numerical  results  are  given  showing  some 
of  the  features  of  the  method. 


1  Introduction 

In  this  paper,  a  novel  three-dimensional  (3D)  hybrid  FEM  and  High  Frequency  (HF)  methodology  for 
the  efficient  analysis  of  general  scattering  problems  is  presented.  It  makes  use  of  FEM  for  the  regions 
with  small  and  complex  features  (arbitrary  geometry  and/or  inhomogeneous  — even  anisotropic —  me¬ 
dia)  and  a  HF  technique  for  the  analysis  of  the  electrically  large  objects  of  the  structure,  taking  into 
account  mutual  interactions  between  the  FEM  domains  and  the  objects  analyzed  with  the  HF  tecnique. 
The  hybrid  methodology  proposed  is  independent  of  the  particular  HF  technique  used.  However,  the  im¬ 
plementation  presented  in  this  paper  corresponds  to  the  use  of  Physical  Optics  (PO)  as  the  HF  technique. 
PO  has  been  chosen  because  of  its  accuracy  when  the  interaction  of  the  FEM  domain  and  the  electrically 
large  objects  is  in  the  near  field  and  because  it  allows  an  easy  implementation  with  distributed  sources 
(as  those  provided  by  FEM). 

The  FEM  part  is  an  improved  implementation  of  the  iterative  methodology  for  general  open-region 
problems  presented  in  [1],  [2],  This  FEM  method  allows  an  easy  hybridization  with  other  techniques. 
Examples  of  hybridization  with  PO  in  the  context  of  scattering  problems  was  presented  in  [3]  and  also  for 
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radiation  problems  in  [4].  In  both  cases,  the  hybridization  was  restricted  to  two-dimensional  problems. 
In  this  paper,  a  3D  implementation  of  this  hybrid  method  is  presented  for  the  first  time.  The  hybrid 
method  is  described  in  Section  2.  Several  numerical  results  of  its  application  are  shown  in  Section  3 
showing  some  of  the  features  of  the  method. 

2  Hybrid  Method 

As  the  hybrid  methodology  proposed  is  independent  of  the  particular  HF  technique  used,  the  description 
of  the  method  made  below  is  general,  i.e.,  valid  for  any  HF  method.  Particularization  to  PO  as  the  HF 
technique  is  made  only  when  necessary. 

Consider  the  problem  illustrated  in  Fig.  1.  It  consists  of  a  region  (bounded  by  an  arbitrarily  shaped 
auxiliary  boundary  S')  with  small  features  and  a  complex  configuration  where  there  are  several  materials 
(which  may  be  anisotropic  and  with  electric  and/or  magnetic  losses)  and  conducting  objects  of  arbitrary 
shape.  A  rigorous  method  such  as  FEM  is  needed  for  the  analysis.  The  region  external  to  S'  consists  of 
an  open  homogeneous  medium  with  the  presence  of  some  electrically  large  objects  (typically  conducting 
and  dielectric  coated  conducting  objects).  The  electrical  size  of  those  objects  makes  appropriate  the  use 
of  an  (approximate)  asymptotic  technique  for  the  analysis  of  that  region.  Thus,  these  objects  will  be 
referred  as  HFBs  (HF  Bodies).  An  incident  field  is  considered  as  excitation. 

The  FEM  domain  (f2)  is  not  truncated  at  S'  but  at  a  boundary  S,  located  outside  S'  and  in  such  a  way 
that  the  region  between  S  and  S'  is  homogeneous.  Boundary  S  may  be  arbitrarily  shaped  but  typically 
it  is  selected  conformal  to  S'.  Distance  from  S'  to  S  is  usually  small,  typically  in  the  range  of  0.05A  to 
0.2A.  Thus,  the  FEM  domain  can  be  truncated  very  close  to  the  sources  of  the  problem. 

The  FEM  formulation  is  based  on  the  vector  wave  equation, 

V  x  (./T1  V  x  V)  -  k20grV  =  0,  (1) 
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V 

fr 

fr 

PEC 

PMC 

E  Formulation 

E 

£r 

Hr 

n  x  V  =  0 

n  x  ^/“1V  xk)=0 

H  Formulation 

H 

flr 

£r 

n  x  (/r_1V  xk)  =0 

n  x  V  =  0 

Table  1:  Correspondences  for  FEM  formulation 


where  k0  is  the  free  space  number  wave,  V  refers  to  E  or  H,  and  fr,  gr  are  defined  accordingly  (see 
Table  1).  The  boundary  conditions  inside  S  may  be  of  Perfect  Electric  Conductor  (PEC)  and  Perfect 
Magnetic  Conductor  (PMC)  that  corresponds  to  Dirichlet  or  Neumann  type  depending  on  the  formula¬ 
tion  used  (see  Table  1). 

A  Cauchy  boundary  condition  is  used  on  S 

n  x  (v  xV^j  +  j  k0n  x  (n  x  v'j  —  ^  (2) 

Differential  equation  (1)  is  discretized  inside  Q  and  a  sparse  system  of  equations  is  obtained 

Kn  KJS 
_  KSi  KSs 

where  the  sub-indexes  S  and  /  refer  to  the  degrees  of  freedom  g  on  S  and  in  the  interior  region,  respec¬ 
tively.  Vector  by  is  a  function  of  the  value  of  the  boundary  condition  on  S,  i.e.,  T  (see  (2)). 

In  order  for  (3)  to  be  solved  by,  i.e.,  T,  must  be  known.  An  initial  value  of  '  is  calculated  as 

#0)  =  fi  x  (V  x  V(0))  +  j  k0n  x  (n  x  V (°})  (4) 


(3) 


—*  (0)  ~*(0)  -*inc  -*inc 

where  V  corresponds  to  the  incident  field  in  the  presence  of  the  HFBs,  i.e.,  V  =  V  +  V HFT  being 

-* inc  ->  inc 

V  the  incident  wave  of  Fig.  1  and  V HFT,  the  incident  wave  field  scattered  by  the  HFBs  — in  absence 
of  FEM  domain — ).  Then  (3)  is  solved.  From  the  FEM  field  solution  and  making  use  of  the  Equivalence 
Principle,  the  value  of  the  scattered  field  (and  its  curl)  on  S  are  computed  from  the  values  of  the  field  and 
its  curl  on  the  auxiliary  boundary  S'.  The  contributions  of  HFBs  are  calculated  using  a  HF  technique, 
i.e., 


VSC(r)  =  j>  -  (V  x  V(f '))  x  VG+  (V  x  (V  x  V(fr 


g+4vvg) 


ds'  +  V 


HFT 


(5) 


V  x  V'c  (f)  =  j  -  (h'  x  (V  x  V  (r 


x  VG  +  (  n  x  V 


/ 


(6) 


where  re  S,  f  '  €  S',  and  G  corresponds  (typically  but  not  necessarily)  to  the  free  space  Green’s 

— *  sc  — *  SC 

function.  Specifically,  in  the  implementation  presented  in  this  paper  Hhft  and  V  x  VrHFT  correspond  to 

->sc  ->  (0)  -*  sc 

the  scattered  field  and  its  curl  computed  by  PO.  The  total  field  is  obtained  as  V  =  V  +  V  and  a 

-  -(i) 

new  value  of 'T,  'k  ,  is  obtained  using  (2). 

In  order  to  compute  VHFT  (and  V  x  Hhft),  the  radiated  field  on  HFBs  from  the  equivalent  sources  on 
S'  is  calculated,  mutual  interactions  between  HFBs  are  applied,  and  the  reradiated  field  (from  HFBs  to 
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boundary  S)  is  computed.  This  process,  although  implemented  as  above,  is  conceptually  equivalent  to 
the  use  of  a  modified  Green’s  function  Grn  (r,  r ' )  to  compute  the  values  of  the  scattered  field  V  (and 
its  curl)  on  S  from  the  equivalent  currents  on  S',  i.e.,  where  Gm  (r,  f ')  =  G  (r,  r  ')  +  GHFT  (r,  r ').  See 
Fig.  2. 

^(i) 

Thus,  the  HF  method  is  easily  hybridized  with  FEM.  The  new  value  is  compared  with  the 

— (0)  -(1)  -*(o) 

previous  one  T  .  In  general,  T  will  be  different  from  T  .  Then,  the  FEM  system  (3)  is  solved  using 
T !l  as  the  value  of  the  boundary  condition  .  The  iterative  process  continues  until  a  certain  error  criteria 
is  satisfied;  typically,  the  7.2- norm  of  the  difference  between  the  values  of  T  between  two  consecutive 
iterations  being  lower  than  a  given  value. 

Fig.  2  shows  the  flow  chart  of  this  iterative  hybrid  technique.  The  method  described  above  has 
several  advantages:  the  FEM  domain  may  be  truncated  with  a  minimum  number  of  unknowns,  the 
original  sparse  and  banded  structure  of  the  FEM  matrices  is  retained,  existing  code  for  non-open  region 
problems  may  be  reutilized,  and  an  easy  implementation  of  the  hybridization  by  means  of  a  simple 
modification  of  the  iterative  FEM  loop  is  obtained.  All  that  is  achieved  at  the  expense  of  performing  a 
number  of  iterations  where  the  numerical  cost  of  the  second  and  subsequent  iterations  is  very  small.  If 
a  direct  solver  is  used  to  solve  the  FEM  problem  at  each  iteration  cycle,  the  factorization  of  the  FEM 
matrix  must  be  performed  only  once  at  the  first  iteration.  Thus,  the  FEM  solution  for  the  second  and 
subsequent  iterations  are  obtained  by  simple  backward  substitution.  If  an  iterative  solver  is  used,  the 
solution  of  the  previous  iteration  cycle  may  be  used  as  an  initial  guess  for  the  next  iteration  of  the  solver. 

It  is  worth  noting  that  the  problem  of  Fig.  1  depicts  a  simplified  example  used  to  describe  the  hybrid 
approach.  In  a  general  case,  several  FEM  domains  may  exist.  Also,  infinite  ground  planes  may  be  taken 
into  account  analytically  using  the  free  space  Green’s  function  with  image  theory.  Numerical  results  are 
shown  in  the  next  section. 

3  Results 

An  example  of  a  dielectric  (er  =  2—j,  //,.  =  1)  coated  metallic  cube  above  a  perfect  conductor  plane  (see 
Fig.  3)  has  been  chosen  to  illustrate  some  of  the  features  of  the  hybrid  method.  The  dielectric  boundary 
is  used  as  the  auxiliary  surface  S'.  The  FEM  domain  boundary  S  is  chosen  conformal  to  S'.  The  FEM 
domain  is  discretized  with  an  irregular  mesh  of  1053  tetrahedra  that  corresponds  to  7694  degrees  of 
freedom  making  use  of  second-order  curl-conforming  tetrahedral  elements  as  those  of  [2],  [5].  For  the 
solution  of  the  FEM  system  a  frontal  solver  from  [6]  is  used.  The  preprocess  (geometry  definition,  mesh 
generation,  etc)  and  postprocess  (field  visualization,  etc)  have  been  made  with  GiD  [7]. 

The  plane  is  taken  into  account  in  the  analysis  by  means  of  PO  where  PO  current  density  JPO  is 
obtained  from  the  incident  magnetic  field  on  the  plane.  The  PO  plane  is  discretized  with  regular  quadri¬ 
lateral  elements  of  side  equal  to  0. 1  A.  Inside  each  element  the  current  is  approximated  by  linear  functions 
which  preserve  its  normal  continuity  between  adjacent  elements.  Fig.  4  shows  \Jpb  \  for  the  first  six  iter¬ 
ations  of  the  hybrid  method  when  a  9  polarized  plane  wave  is  illuminating  the  plane  (of  W  =  2  A)  from  a 
normal  direction.  It  may  be  observed  how  the  current  is  modified.  Initial  current,  \JPq\,  is  constant  and 
its  value  is  2/r/0.  fpQ 1 1  —  JpQ |  diminishes  with  each  iteration  and  \Jpo\  converges  toward  a  function 
almost  constant  on  the  illuminated  surface  with  small  amplitude  in  the  shade  region.  Furthermore,  a 
small  orthogonal  polarization  component  is  generated.  Fig.  5  shows  the  bistatic  RCS  results  with  W 
as  parameter.  As  the  total  scattered  field  increases  with  W,  the  scattered  field  shown  in  the  figure  is 
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the  total  scattering  fields  minus  the  reflected  field  by  the  plane  when  there  is  no  FEM  object  (the  cube). 
The  results  for  the  latter  are  obtained  by  the  iterative  FEM  solution  using  dyadic  Green’s  function  with 
Dirichlet  (Neumann)  condition  for  E-formulation  (H-formulation).  It  may  be  observed  how  the  results 
obtained  with  the  hybrid  method  for  the  finite  size  plane  converge,  as  the  plane  width  W  is  getting  larger, 
to  those  obtained  with  FEM  and  and  an  infinite  width  plane.  The  results  mentioned  above  have  been 
obtained  using  the  E-formulation.  Analogous  results  have  been  obtained  with  the  H-formulation.  The 
number  of  iterations  of  the  method  is  around  25  for  ||6^  —  6^  ^  ||  <  10-4. 

4  Conclusions 

A  hybrid  methodology  for  the  efficient  analysis  of  scattering  problems  where  small  and  complex  features 
are  combined  with  electrically  large  objects  has  been  presented.  It  makes  use  of  FEM  for  the  regions 
with  small  and  complex  features  and  an  HF  technique  for  the  analysis  of  the  electrically  large  objects  of 
the  problem.  Specifically,  its  3D  implementation  with  PO  as  the  HF  tecnique  is  presented  and  numerical 
results  are  shown  to  illustrate  the  main  features  of  the  approach. 
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Fig.  2:  Flow  chart  of  the  iterative  hybrid  FEM-FIF  Technique 
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A 

Fig.  3:  Dielectric  coated  cube  cube  above  a  perfect  conducting  plane 
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e)  Fifth  iteration 
Fig.  4:  |«7po|  on  plane  (IF 


f)  Sixth  iteration 

=  2  A)  versus  iteration  number  for  problem  of  fig.  3 
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Bistatic  Radar  Cross  Section.  elnc=0Q,  olnc=0- 


a)  <^  =  0° 


Bistatic  Radar  Cross  Section.  elnc=0e,  olnc=0- 


h)  (j>  =  90° 

Fig.  5:  Bistatic  RCS  for  problem  of  fig.  3  for  different  W 
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